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Abstract. The psi function 4){x) is defined by ^{x) = and (x) for 

i g N denote the polygamma functions, where r(x) is the gamma function. In 
this paper we prove that a function involving the difference between [■(/>' (x)]^ + 
tp"{x) and a proper fraction of x is completely monotonic on (0, oo). 



1. Introduction 

We recall [3, Chapter XIII] and [4, Chapter IV] that a function / is said to be 
completely monotonic on an interval / if / has derivatives of all orders on / and 



for X E I and n > 0. 

The famous Bernstein- Widder Theorem (see [4, p. 160, Theorem 12a]) states 
that a function f{x) on [0, oo) is completely monotonic if and only if there exists a 
bounded and non-decreasing function a{t) such that 



Jo 

converges for x G [0, oo). This says that a completely monotonic function f{x) on 
[0, oo) is a Laplace transform of the measure a{t). 

We also recall that the classical Euler gamma function r(a;) is defined by 
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The logarithmic derivative of T{x), denoted by ^{x) — ^j^^^j , is cahed the psi or 
di-gamma function, and the derivatives V'*'*' (x) for i £ N are respectively called the 
polygamma functions. In particular, the functions ^' {x) and i^" {x) are called the 
tri-gamma and tetra-gamma functions. 

In [2, p. 208, (4.39)], it was established that the inequality 

holds for X > 0, where 

p{x) = 75a;^° + 900a;^ + 4840a;* + 15370a;^ + 31865a;^ + 45050a;^ 

(5) 

+ 44101a;^ + 29700a:^ + 13290x2 + 3600x + 450. 

The aim of this paper is to prove the complete monotonicity of the difference 
between two functions on both sides of the inequality (4). 
Our main result may be stated as the following theorem. 

Theorem 1. The junction 

,(x) = [/(xr + /'(x)-^^^^^^ (6) 

is completely monotonic on (0, oo), where the function p{x) is defined by (5). 

Remark 1. By the definition of completely monotonic functions, it is easy to see 
that Theorem 1 is stronger than the inequality (4), so Theorem 1 generalizes the 
inequality (4). 

2. Proof of Theorem 1 

By the recursion formula [1, pp. 258 and 260, 6.3.5 and 6.4.6] 

(.T + 1) = {x) + (7) 

for X > and n € N, direct calculation produces 

g{x) - g{x + 1) = [V'(x) - V'(ar + 1)] [^'(x) + V'(ar + 1)] 

p{x) p{x + 1) 



+ [V'"(x)-V"(x + l)] 



900x4(x + l)io 900(x + l)4(x + 2)io 



1 ■ 


2 


p{x) 


p{x + 1) 


x^ 


x^ 


900x4(x + 1)10 


900(x+ l)4(x + 2)io 
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2a;2 X 1800a;2(a; + l)io 1800(a; + l)4(a; + 2)1° 

■ , 1 3 28 51 251 



2x 4a;2 5(a; + 1) 10(a; + 2) 120(a;+l)2 



X 

331 7 17 13 49 



120(a; + 2)2 6(a; + l)3 12(a; + 2)3 90(x+l)4 72(a; + 2)4 
13 47 1 1 



+ 



180(a:+l)5 180(x + 2)5 120(a;+l)6 60(a; + 2)6 
1 2 1 13 



180(a; + l)7 45(a; + 2)7 200(a; + 1)8 600(a; + 2)8 
1111 

7^7^777 ! TTn "7^777 ! TTTn T7777777 ! TTTn" 



900(a;+l)9 450(a; + 2)9 1800(a;+l)io 450(x + 2)io 



Using the formula [1, p. 255, 6.1.1] 

1 1 



j.r— 1 xt 



dt 



r(r) Jo 

for r > and a; > and the integral representations [1, p. 260, 6.4.1] 

poo j.n 

^(")(a;) = (-1)"+M -e-^^dt 

Jo 1 - e * 

for n G N and x > gives 

../X ( i 1 3 28 , 51 2t 251 , 

ma;) = / 1 e~* H e"2* H te"* 

io Vl-e-* 2 4 5 10 120 

331 oj. 7 o + 17q o+ 13q j. 49q f>j. 

+ ie"^* re"* H re"^* H ^e"* H fe~^* 

120 12 24 540 432 



4320 4320 14400 7200 

+ I +6 -t ?_+6 -2t , 5^ .7 -L ]^ .7 -2* 

129600 32400 1008000 3024000 

+ I t^e-* - t^e-^* - fe-* 

36288000 18144000 653184000 

1 



163296000 



1 1 

■- / [I63296000(i - 2)e^* - e^* (t^ - 18t^ 

653184000 7o - 1 

- 648i^ - 5040i^ + 45360f^ + 1965600f'* - 15724800t^ 

+ 381024000*^ - 1856131200* + 3331238400) + e*(5i^ - 54f8 

- 3456i^ - 45360t'' - 45360t^ + 9072000t^ + 58363200*^ 
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0"{t) 



+ 843696000t2 + 435456000t+ 6989068800) - 4(t^ - 9f 
- 702f - lOOSOt*^ - 22680i^ + 1776600f'^ + 18522000^^ 
+ 115668000*2 + 450424800* + 832809600)] e-(^+2)* d* 



163296000(3* - 5)e^* - (4806345600 - 2950214400* 

+ 714873600*2 - 23587200*''' + 4158000** + 60480*-" - 14616**^ 

- 1440*^ - 27** + 2*^)e2* + (7424524800 + 2122848000* 

+ 1018785600*2 + 94651200*^ + 8845200*"^ - 317520*^ - 69552**^ 

- 3888*^ - 9** + 5*^) e* - 36(50047200 + 25704000* + 6174000*^ 
+ 789600*3 _ 12600*"* - 6720*^ - 546*^ - 8*^ + **), 
489888000(3* - 4)6^* - 4(1665619200 - 1117670400* 

+ 339746400*2 - 7635600*^ + 2154600** + 8316*^ - 9828*^ 

- 774*^ - 9** + *^)e2* + (9547372800 + 4160419200* 
+ 1302739200*2 + 130032000*^ + 7257600** - 734832*^ 

- 96768**5 _ 3950^7 ^ 35^8 ^ 5^9^ _ 144(6426000 + 3087000* 
+ 592200*2 - 12600*^ - 8400** - 819*^ - 14**5 ^ 3*^) , 
4408992000(*- l)e^* - 4(2213568000- 1555848000* 

+ 656586000*2 _ 6652800*^ + 4350780** - 42336*^ - 25074*'' 

- 1620*^ - 9** + 2*9) e2* + (13707792000 + 6765897600* 
+ 1692835200*2 + 159062400*^ + 3583440** - 1315440*^ 

- 124488*^ - 3672*^ + 81** + 5t^)e* - 1008(441000 + 169200* 

- 5400*2 - 4800*^ - 585** - 12*^ + 2*^) , 
4408992000(3* - 2)6^* - 16(717822000 - 449631000* 

+ 323303400*2 + 1024380*^ + 2122470** - 58779*^ - 15372**5 

- 828*^ + *9)e2* + (20473689600+ 10151568000*+ 2170022400*2 




Straightforward computation yields 
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+ 173396160t^ - 2993760t^ - 2062368t^ - 150192^^ - 3024i'' 
+ 126t^ + 5t^)e* - 12096(14100 - 900i - UOOt^ - 195t^ - Si"^ + 1^), 
e^^\t) = 13226976000(34- l)e^* - 16(986013000 - 252655200i 
+ 649679940*^ + 10538640*^ + 39510454^ - 209790i^ 

- 36540t^ - 1656t^ + 9t^ + 2t^)e^^ + (30625257600 

+ 144916128004+26902108804^ + 1614211204^ - 133056004'' 

- 29635204^ - 1713604*' - 20164^ + 1714^ + 54^) e* 

- 60480(-180 - 4804 - 1174^ - 44^ + 4*), 
6l(6)(t) = 11904278400046^* - 64(429842700+ 1985123704 

+ 3327439504^ + 92203654^ + 17132854"^ - 1597054^ 

- 211684*^ - 8104^ + 94^ + 4^)6^* + (45116870400 

+ 198720345604 + 31744742404^ + 1081987204^ - 281232004* 

- 39916804'^ - 1854724^' - 6484^ + 2164^ + 54^^)6' 

- 120960(-240 - 1174 - 64^ + 24^), 

6»(^'(4) = 119042784000(1 + 34)e^* -64(1058197770+ 10625126404 

+ 6931489954^ + 252938704^ + 26280454* - 4464184^ - 480064^ 

- 15484^ + 274^ + 24'^) e^' + (64988904960 + 262209830404 
+ 34990704004^ - 42940804''' - 480816004* - 51045124^ 

- 1900084^ + 10804^ + 2614^ + 54^) e* - 362880 (24^ - 39 - 44), 
6»(8)(4) = 357128352000(2 + 34)e^* - 128(1589454090+ 17556616354 

+ 7310898004^ + 305499604^ + 15120004* - 5904364^ - 534244^ 

- 14404^ + 364^ + 24^) e^* + (91209888000 + 332191238404 
+ 34861881604^ - 1966204804^ - 736041604* - 62445604^ 

- 1824484^ + 31684^ + 3064^ + 54^)e* - 1451520(4 - 1), 
e'^^\t) = 3214155168000(1 + 4)e^* - 128(4934569815 + 49735028704 

+ 15538294804^ + 671479204^ + 718204* - 15014164^ - 1169284^ 
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- 2592f^ + 90<® + 4t9)e2* + (124429011840 + 40191500160t 
+ 2896326720i2 - 491037120f3 - 104826960f'' - 7339248f^ 

- 160272i*^ + 5616i^ + 351t^ + 5t^)e* - 1451520, 

6»(io)(t) = e* [164620512000+ 45984153600* + 1423215360i^ - 910344960*^ 

- 141523200*'* - 8300880t^ - 120960*'^ + 8424f^ + 396t** + 

- 512(3710660625 + 3263666175^ + 827275680i2 + 33645780i3 

- 1840860*^ - 926100f^ - 63000f*' - 1116t^ + 54i^ + 2t^)e* 
+ 3214155168000(4+ 3i)e2*] 

0[{t) = 3214155168000(11 + 6t)e^* - 512(6974326800 + 4918217535f 

+ 928213020*^ + 26282340*^ - 6471360t^ - 1304100*^ - 70812i^ 

- 684i^ + 72t^ + 2t^)e* + 9(5109350400+ 316270080t 

- 303448320f2 - 62899200i3 - 4611600*'* - 80640i^ 
+ 6552**^ + 352*^ + 5**), 

e'l{t) = 8 [1607077584000(7 + 3*)e2* - 64(11892544335 + 6774643575* 

+ 1007060040*2 + 396900*^ - 12991860*^ - 1728972*^ - 75600*^ 

- 108*^ + 90** + 2*^)e' + 9(39533760 - 75862080* - 23587200*^ 

- 2305800*3 - 50400*^* + 4914*^ + 308**^ + 5*^)] , 

e^^\t) = 8 [1607077584000(17 + 6*)e2* - 64(18667187910 + 8788763655* 
+ 1008250740*2 - 51570540*^ - 21636720*'' - 2182572*^ 

- 76356*^ + 612*^ + 108** + 2t^)e* - 63(10837440 + 6739200* 
+ 988200*2 + 28800*3 - 3510*^ - 264*^ - 5**^)] , 

e[^\t) = 16 [3214155168000(10 +3*)e2* - 32(27455951565 

+ 10805265135*+ 853539120*2 - 138117420*^ - 32549580*"* 

- 2640708*^ - 72072**^ + 1476*^ + 126** + 2*^)e* 

- 945(224640 + 65880* + 2880*2 - 468*^ - 44*"* - t^)] , 
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0[^\t) = 16 [3214155168000(23 + 6i)e2* - 32(38261216700 

+ 12512343375* + 439186860*2 - 268315740*^ - 45753120*^ 

- 3073140*^ - 61740*'^ + 2484*^ + 144*^ + 2t^)e* 

- 945(65880 + 5760t - U04t^ - 176t^ - 5t^)] , 
0^^\t) = 64 [3214155168000(13 + 3t)e2* - 8(50773560075 

+ 13390717095t- 365760360t2 - 451328220t3 - 61118820*4 

- 3443580*^ - 44352i<^ + 3636t^ + 162t^ + 2f)e* 

- 945(1440 - 702t - 1321"^ - 5t^)] , 

e^^\t) = 64[3214155168000(29+ 6<)e2* - 8(64164277170 

+ 12659196375t- 1719745020^^ - 695803500t3 - 78336720t'* 

- 3709692*^ - 18900**^ + 4932t^ + 180t^ + 2f)e* 
+ 2835(234 + 88t + 5^2)], 

e^^\t) = 128 [6428310336000(16 + 3f)e2' - 4(76823473545 

+ 9219706335^ - 3807155520*^ - 1009150380*^ - 96885180*"* 

- 3823092*^ + 15624**^ + 6372*^ + 198*^ + 2f)e* + 2835(44 + 5*)] , 
e^^\t) = 128 [6428310336000(35 + 6*)e2* - 4(86043179880 

+ 1605395295* - 6834606660*^ - 1396691100*^ - 116000640*^ 

- 3729348*^ + 60228*^ + 7956*^ + 216*^ + 2f)e^ + 14175] , 
e[^°\t) = 512e* [6428310336000(19+ 3*)e* - 87648575175 + 12063818025* 

+ 11024679960*2 + 1860693660*^ + 134647380*^ + 3367980*^ 

- 115920*^ - 9684*^ - 234** - 2*^] 
= 512e*02(*), 

6»2(*) = 9 [1340424225 + 2449928880*+ 620231220*^ + 59843280*^ 

+ 1871100*4 - 77280*^ - 7532*^ - 208*^ - 2** 

+ 714256704000(22 + 3*)e*] , 
e'^it) = 72 [306241110+ 155057805*+ 22441230*2 + 935550*^ 
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- 48300t^ - 5649t^ - 182t^ - 2f + 89282088000(25 + 3t)e*] , 
e^^\t) = 504[22151115 + 6411780* + 400950*^ - 27600*^ - 4035*^ 

- 156t^ - 2t^ + 12754584000(28 + 3t)e*] , 
e^^\t) = 6048 [534315 + 66825t - 6900*^ - 1345f^ - 65t^ - 

+ 1062882000(31 + 3t)e*], 
ei^\t) = 30240 [13365 - 2760t - 807t^ - 52t^ - 

+ 212576400(34 +3t)e*], 
ef\t) = 60480 [106288200(37 + 3t)e* - 1380 - 807t - TSt^ - 2t^] , 

(t) = 181440 [35429400(40 + 3t)e' - 269 - 52t - 2f] , 
0^^\t) = 725760 [8857350(43 + 3t)e* - 13 - , 
e^^\t) = 725760 [8857350(46 + 3i)e* - l] . 

is easy to calculate that 



^'(0) 


= 0, 


^"(0) 


= 0, 


^(3)(0) 


= 0, 


^W(O) 


= 0, 


^(^)(o) 


= 1632960000, 


0(*^)(O) 


= 17635968000, 




= 116321184000, 


0(8) (0) 


= 602017920000, 


e(9)(0) 


= 2706957792000, 


0(1°) (0) 


= 11121382944000, 


e[{0) 


= 31830835680000, 


0'I{O) 


= 83910208435200, 




= 208999489144320, 




= 500203983121920, 




= 1163218362768000, 


0f\Q) 


= 2648180949926400, 




= 5932619924353920, 


0f\o) 


= 13125845965639680, 




= 28754776198995840, 




= 62489726878118400, 




= 141434891210025, 


om 


= 160729807759920, 




= 180003853569960, 




= 199280851953120, 


(0) 


= 218562955581600, 




= 237847398969600, 




= 257132364632640, 


0i'\o) 


= 276417335013120, 
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6»(^) (0) = 295702274730240. 

Since 0^\t) is increasing, so 0^\t) > on (0, oo), which means that 0^\t) is 
increasing and positive on (0, oo). By the same argument, it is derived that the 
functions e^\t) for 1 < i < 9, e^\t) and e^'\t) for 1 < i < 10 are mcreasing 
and positive on (0,oo). Therefore, the function 6{t) is increasing and positive on 
(0,oo), which imphes that the function H{x) is completely monotonic on (0, oo). 
Because the function ^ is completely monotonic on (0, oo) and the product of finite 
completely monotonic functions are also completely monotonic, we obtain that the 
function g{x) — g{x + 1) is completely monotonic on (0, oo), which is equivalent to 

{-l)''[g{x) - g{x + 1)]W = (-l)'=5W(a;) - (-l)'=sW(x + 1) > 

for < on (0, oo). By induction, we have 

(-1) V^a;) > (-l)'=ffW(a; + 1) > (-1) V^a^ + 2) > ■ • ■ 

> {-l)''g^''\x + m) > • • ■ > lim [{-1)'' g^''\x + m)] = 

for /c > on (0, oo). The proof of Theorem 1 is complete. 

Remark 2. Simplification yields that the function 

"M=^'W- 1800.'(l+lr(2 + .f° 

is completely monotonic on (0, oo), where 

Q{x) = 1382400 + 21657600a; + 162792960a;2 + 778137600a;^ 
+ 2645782983a;^ + 6789381590a;^ + 13626443025x^ 
+ 21889330810x^ + 28579049475a;^ + 30634381522.t'' 
+ 27125436630x1° + 19896883200x" + 12088287630a;i2 
+ 6063596590a;i^ + 2494770300a;i* + 8329584000;^^ 
+ 222060150x1^ + 46134540x1^ + 7195500x1^ 
+ 792300x1^ + 54900x2° + 1800x^1 



for X G (0, oo). 
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